Dynamical multistability in high-finesse micromechanical optical cavities 
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We analyze the nonlinear dynamics of a high-finesse optical cavity in which one mirror is mounted 
on a flexible mechanical element. We find that this system is governed by an array of dynamical 
attractors, which arise from phase-locking between the mechanical oscillations of the mirror and the 
ringing of the light intensity in the cavity. We describe an analytical approximation to map out the 
diagram of attractors in parameter space, derive the slow amplitude dynamics of the system, includ- 
ing thermally activated hopping between different attractors, and suggest a scheme for exploiting 
the dynamical multistability in the measurement of small displacements. 
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Introduction. - The radiation pressure exerted by light 
stored in an optical cavity couples the mechanical de- 
grees of freedom of the cavity mirrors to the optical field. 
The dynamics which result from this coupling offer novel 
means for manipulating both the light and the mirrors 
themselves, for example by generating squeezed [l| or 
entangled states of light |2| or matter |3| , or by tailoring 
the mechanical properties of the mirrors 0, IM H • 
The radiation pressure also sets fundamental and prac- 
tical limits on the sensitivity of a range of experiments, 
from the km-scale optical cavities in gravitational wave 
observatories such as LIGO to the /im-scale cavities used 
to probe the motion of micromechanical systems, e.g. by 
enforcing the standard quantum limit for displacement 
measurements p| llfll ITU llij or introducing mechanical 
instabilities [lElHuB- 

To date the theoretical description of these systems 
has been almost entirely in terms of linearized equations 
of motion, valid for small mirror oscillations. While this 
approach has been adequate to predict a wide range of 
novel effects, it neglects the inherent nonlinear ity of these 
devices and cannot describe their behavior in dynami- 
call y u nstable regimes. As advances in micromachining 
EES 111 Gil and optical cavity fabrication (UHH pro- 
duce devices with stronger optomechanical coupling and 
weaker damping than was previously possible, the non- 
linearities and dynamical instabilities of these systems 
will become increasingly important. Here we present an 
analytic theory which provides a full description of the 
nonlinear regime. We find that the dynamics is domi- 
nated by an array of stable dynamical attractors which 
correspond to phase-locking between the mirror motion 
and the ringing in the cavity optical field. These re- 
sults are confirmed in detail by numerical simulations of 
the microscopic equations of motion. The theory pre- 
sented here also describes the statistical behavior of the 
cantilever oscillation amplitude and thermally-activated 
transitions between the dynamical attractors. 

The model. - We consider the setup of Fig. The 
coupling between the cantilever position and the light in- 
tensity of a given cavity mode is described by the Hamil- 
tonian H Cdcv = Hul(1 — x/l)a^a, provided we may neglect 
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Figure 1: The setup analyzed in the text (not drawn to scale). 



the finite travel-time of the light between the mirrors. 
The cavity of length / is resonant with the incoming laser 
radiation of frequency uol when x = 0, and we assume the 
excursions in x to be small enough to avoid the other res- 
onances. This coupling gives rise to a radiation pressure 
force F = (hjL/l)o)a. In the following we will consider 
the purely classical nonlinear dynamics, where a is re- 
placed by the coherent light amplitude a. We use the 
ring-down time of the cavity 7 -1 , the resonance width 
Sx = Ij/ujl: and the cantilever mass m as convenient 
new units, and rescale a such that it becomes 1 at reso- 
nance: a = a OTlg e luJLt / -\/n max , where the maximum pho- 
ton number n max = 4Pi n / (jHcol) is linear in the input 



power Pi n . Then the coupled equations of motion for a 
and x read 



(i) 

(2) 



V\a\ 2 — Wq(x — xq) — Tx 



The oscillator frequency coo and mechanical damping rate 
T are fixed, while the detuning from resonance xq may be 
controlled, either by changing the laser frequency or ap- 
plying a static force to the cantilever. All the other con- 
stants are combined into the dimensionless input power 



V 
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m^H 2 mSx 2/ y 
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We emphasize that the scale of the radiation pressure, 
set by V, grows with the fourth power of the cavity fi- 
nesse (V cx 7 -4 with 7 = (c/2/)T, where T is the mirror 
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transmission). The nonlinear effects we investigate be- 
come important for T <C V . 

This work is motivated in part by the possibility of 
using the fabrication process described in [H Eg to 
integrate ultrasensitive cantilevers with high-reflectivity 
dielectric mirrors, and to use these mirrors as part of 
a high- finesse optical cavity at low temperatures. In 
[HI, cantilevers with spring constants of ~ 10 _3 N/m 
and mass 10 _11 kg were integrated with lithographically- 
patterned high-quality heterostructures similar in design 
to high-reflectivity dielectric mirrors. These cantilevers 
were then used in an optical interferometer with an in- 
cident power Pi n ~ 20nW at T ~ 0.3 K. By making 
modest adjustments to the dimensions of these devices, 
it should be possible to realize cantilevers with spring 
constant ~ lN/m, mass m ~ 10 _10 kg, and mechanical 
Q ~ 10 5 supporting a mirror capable of achieving a cavity 
finesse of 7r/T ~ 10 5 , where I ~ 10cm. Assuming these 
values, the cavity ringdown rate 7 ~ 10 5 Hz becomes com- 
parable with the cantilever frequency (c^o ~ 1 in units 
of 7), and the dimensionless damping rate T ~ 10 -5 is 
much smaller than V ~ 1. In contrast, our estimate for 
the recent low-finesse experiment on radiation cooling of 
a cantilever is roughly V ~ 10 -18 <C T (both for 
radiation-pressure and photothermal effects). 

The cavity resonance peak a(x) = 1/(1 — 2ix) gives rise 
to a rounded step-like barrier in the effective static can- 
tilever potential obtained by integrating the right-hand- 
side (rhs) of Eq. 0, 



V^(x) = ^(x- 



x ) 2 



arctan(2x) 



(4) 



Depending on the parameters (c^c^o,^), there are one 
or two local minima of T4ff, opening the possibility of 
static bistability 0,0- However, it is known 4] that the 
time-lag of the radiation pressure force generated by the 
finite cavity ring-down time 7 -1 introduces additional 
damping or anti-damping to the left or right of the bar- 
rier, respectively. Here we will focus on the regime where 
the anti-damping leads to a linear dynamical instability 
discovered previously 0, 0, Q, , such that there is no 
stable stationary solution x = 0. Then the system set- 
tles into self-sustained oscillations, whose full nonlinear 
dynamics we explore here. 

Dynamics in the unstable regime. - Direct numerical 
solution of Eqs. @ reveals that x(t) carries out 

approximately sinusoidal oscillations at the unperturbed 
frequency cjq, 



x(t) 



Acos(uiot) . 



(5) 



This fact is essential in allowing us to derive an analyt- 
ical theory for the nonlinear cantilever dynamics. The 
light amplitude a(t) develops a more complicated be- 
haviour (Fig. |2|). It experiences a sharp rise each time 
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Figure 2: Left: Light intensity oscillations for given sinu- 
soidal cantilever motion (with uo = 1, x = xo = 5, A = 20). 
Right: Resulting dynamical potential U(A), for V = 1 and 
T = 10" 4 , 10" 3 , 10" 2 , 10" 1 (bottom to top curve). Insets 
show |a| 2 (t), with time on the horizontal axis spanning one 
period and the vertical scale equal to 0.4 in each case. Note 
the additional oscillation in \a\ 2 for each subsequent potential 
well, i.e. each new dynamical attractor. 



the cantilever swings through the resonance at x = 0. Its 
dynamics resembles that of a driven damped oscillator 
which is swept through resonance non-adiabatically. The 
exact solution for a given x(t) (Eq. JSJ) can be written 



as a Fourier series, a(t) = e t(p "' a n e 



O inujot 



with 
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and a global phase <p(i) = (A/uq) sin(cjo^), where J 1 
the Bessel function of the first kind. 

Dynamical multist ability. - We can determine the pos- 
sible dynamical attr actors (x, A) of the cantilever motion 
by imposing two conditions that follow from Eq. JJJ for 
any periodic motion. The total time-averaged force (x) 
has to vanish, and the net power input via the radiation 
pressure force must equal the power dissipated through 
friction, (xx) = 0. This translates into 

V{\a(t)\ 2 ) = u 2 (x-x ) (7) 

Prad = V (Ht)\ 2 x) = Pfrfc = r (X 2 ) . (8) 

The exact dependence on x,i and uoq follows directly 
from the coefficients of Eq. ©: 

<M 2 > = 5>„| 2 (9) 

n 

Prad = (\a\ 2 %) = Aj Im^a*a n+ i (10) 

n 

The power balance equation can be recast into the form 



P ia , d (x,A) _ r 



(11) 
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Figure 3: Density plot of the ratio of reduced radiation power 
input and frictional power loss, P Y &d(x(xo, A), A) / Pf T - lc (A) , 
in the (xo, A)-plane. The contour lines indicate dynami- 
cally stable cantilever oscillation amplitudes A, according to 
Eq. (fTT)) . White contour lines display the approximation 



: xq. Contours are drawn for T/V = 10 , 10 , 10" 



,10" 



(with V = ujo = 1). Red dots show the long-time limit 
(tsim = 2 • 10 4 ) of the amplitude A = x ma x — Xmin in the nu- 
merically exact solution of the original equations of motion, 
obtained for the same values of T/V and a set of random 
initial conditions. 



with Pf r i c (^4) = 0OqA 2 /2. Stable attr actors are those 
where the ratio decreases for increasing A. 

After solving the force balance equation J7J) for x = 
x(xq,A), we can plot the contour lines of the lhs of Eq. 
fnj in the (xq, A)-plane. These trace out the diagram 
of possible equilibrium values of the cantilever oscillation 
amplitude A as a function of detuning xq (Fig. |3J). The 
red dots in Fig. 03 show the results of a (much more 
time-consuming) full simulation of the initial equations 
of motion, Eqs. (CP),®, confirming that our analytical 
approximation works extremely well. 

The initial onset of the instability at some threshold 
Xq < (with Xq = — 1 in Fig. [3J) can be described by 
expanding Eq. J7J) for small A (and x « 0), which is valid 
in the limit r«P. This leads to A « [(1 + 4o;g)(l + 
UqXq/V)/2} 1 / 2 . The position of the onset, x% = -V/ufa, 
and the rise of A do not depend on T in this limit. 

Near the onset we have \x\ <Cxo. However, for the rest 
of parameter space (#o, A), we find that x ~ xq is a good 
approximation (this approximation is shown as the white 
lines in Fig. 03 as compared to the black lines), becom- 
ing exact for V — > 0. Within this approximation, which 
we will adopt in the following, we only need to fulfill the 



power balance equation JSJ). Its lhs does not depend on 
the damping T or the dimensionless light power V, and 
therefore the equilibrium values of A, at a given xo, de- 
pend only on the ratio T/V within this approximation. 
At large A/ouq, we find 
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^4 
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(12) 

This reveals explicitly the periodic structure observed in 
Fig. 03 with periods of ttcoq and 2coq in A and xq, respec- 
tively. Physically, an increase in A enlarges the time- 
dependent detuning x(t) which is equal to the frequency 
of the ringing in the light intensity |«| 2 (t) (Fig. |2J). The 
appearance of the array of attractors is thus due to a 
phase-locking phenomenon, where the oscillations of the 
light intensity (depending on A) seek to be commensurate 
with the fundamental cantilever period. For the exper- 
imental parameters mentioned above, the periods in A 
and xo will be on the picometer scale. 

Amplitude equation of motion. - We can determine the 
slow dynamics of the cantilever oscillation amplitude A 
by equating the change in total cantilever energy E = 

^-A 2 to the net power input: 



dE 
~dt 



— -ftotal — Pr 



rad 



fric 



(13) 



This can be used to obtain the overdamped motion of A, 



dA 



-r- = ^rPt 



Awl 



total 



(14) 



where we have introduced the effective potential for A, 



U(A) 



A 



2V 
~T~ 



f 



A dA' 



(15) 



The first part of U(A) is the oscillator potential itself, 
which drives A towards zero in the absence of radiation 
pressure (V = 0). The second part produces a decaying 
oscillating component, as follows from Eq. (|12|). The 
total potential U(A), obtained from P ra d(A), in general 
displays several local minima, corresponding to the dy- 
namical attractors, before the quadratic rise of the first 
term takes over at large A (see Fig. El right). Solving 
Eq. (|14|) and comparing to ab-initio simulations of the 
full dynamics yields a good agreement, which can be im- 
proved at small \xq\ by using x = x{xq, A) from Eq. 
instead of x = xq. This was also used for Fig. 21 

Stochastic dynamics and Boltzmann distribution. - We 
can model the stochastic dynamics of this driven non- 
linear system within our amplitude equation approach. 
The heat bath responsible for the intrinsic damping of 
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Figure 4: Boltzmann distribution w(A) of the cantilever os- 
cillation amplitude A, as a function of detuning xo, for a 
reduced temperature T = 10. Other parameters: V — cjo = 
1, r = 10 -3 . Note the sharp drop in w(A) at A = —xo 
(dashed blue line), the narrowing (cooling) of the distribu- 
tion for xo < 0, the onset of instability at Xq, and the transi- 
tion (s) above xo = 20, back to the thermal distribution of an 
unperturbed oscillator. Inset: Underlying effective potential 
U (A) (top truncated) , with bare oscillator parabolic potential 
ooqA 2 /2 indicated in red. 



the cantilever also produces a fluctuating force, which 
can be added to the rhs of Eq. J2J) as a white-noise 
term £x(t). The fluctuation-dissipation-theorem de- 
mands (£r(£)6r(0)) = 2TT5(t), with a reduced temper- 
ature T = kBT/(mSx 2 j 2 ). This gives rise to a fluctu- 
ating power £e — x£,x in Eq. (|13|). which in turn leads 
to a stochastic force £a(£) = £e (£)/ C^o ) on tne rns °f 
the equation of motion ([Tl)) for A, with (^(t)^(0)) = 
(TT/ujq) S(t). The cantilever will thus settle into a dy- 
namical, driven equilibrium described by the Boltzmann 
distribution (Fig. HJ) 



w(A) oc Aexp 



U(A) 



(16) 



The prefactor A accounts for the phase space associated 
with the two quadratures of motion. Formally, it is pro- 
duced by an additional deterministic term TT/(2uoqA) 
on the rhs of Eq. ((14)). which stems from rewriting the 
Langevin-version of (p"3|) in Ito-form. Comparison with 
a Langevin-Runge-Kutta [2^ simulation of the original 
Langevin equations for x and a shows good agreement. 
Shot noise may be neglected for the high photon numbers 
of interest here (n max ~ 10 6 ). 

As the scale AU of the corrugations in the oscil- 
latory "washboard-potential" U(A) is set by the ratio 
V/Tj see Eq. (T^J), it can become much larger than the 
light-induced contribution to the static potential V e f[, 
Eq. Q. which scales like V. Thus, at small intrinsic 
damping T, the dynamically unstable cantilever may dis- 
play a strong confinement of the oscillation amplitude 
A in the wells of the dynamical potential, even while 
radiation effects on the static Boltzmann distribution 



exp[— V G fi(x)/T] would be negligible. In our numerical 
example, the scale AU corresponds to a real temperature 
in the tens of Kelvin, below which the Boltzmann distri- 
bution clearly shows the structure of multiple dynamical 
attractors. The Kramers rate for thermal hopping be- 
tween wells is r" 1 = (TAU/2ttu;^) exp[-Af7/T]. Both at 
xo < and large xo > the global minimum of U(A) 
is at A = 0, and w(A) returns to the distribution of the 
unperturbed oscillator. 

Photothermal (bolometric) forces can be another 
source of coupling between light field and cantilever and 
were used for cooling a cantilever 0, Q • Due to the finite 
heat conductance of the cantilever, they feature an ad- 
ditional time-lag r between light intensity and force, in 
addition to the time-lag between cantilever motion and 
intensity already present in the case of radiation pressure. 
Such a force enters the rhs of the cantilever equation of 
motion Q in the form V't' 1 dt f \a\ 2 (t f ) exp[-(t - 
t')/r\, where we have considered the special case r = 
up to now. The prefactor V' is equal to V multiplied by 
the ratio between bolometric and radiation pressure force 
at constant light intensity. The whole preceding analysis 
still applies, provided we divide the sum in Eq. (|TU|) by 
the factor 1 + iujor and add up all contributions to P ra d 
and U(A). The oscillatory structure of the potential is 
fully retained. However, for large ujqt, the total magni- 
tude of this bolometric contribution is suppressed, such 
that only the radiation pressure effects remain. 

Hysteresis and 'latching' measurements. - The pres- 
ence of multiple locally stable dynamical attractors nat- 
urally leads to hysteresis upon sweeping an external pa- 
rameter. In addition, the sensitivity to parameter per- 
turbations is greatly enhanced near transitions between 
attractors. This could lend itself to the convenient mea- 
surement of very small displacements (i.e. small pertur- 
bations 5xo, on the sub-picometer scale): Sweeping back 
and forth in xo can leave the cantilever in either of two 
attractors, depending on whether Sxo was large enough. 
Afterwards, the resulting (stable) amplitude may be mea- 
sured. A similar 'latching' scheme has been proposed and 
implemented in a very different context, using the dy- 
namical bist ability of a strongly driven Josephson junc- 
tion 0. 

Conclusions. - We have studied the nonlinear classical 
dynamics of high- finesse micromechanical optical cavities 
and discovered a form of dynamical multistability that 
can arise in high-finesse cavities with a large mechanical 
quality factor of the oscillating mirror. The cantilever un- 
dergoes self-sustained oscillations whose amplitude set- 
tles into one of several attractors, which we have mapped 
out in parameter space. We have derived an analytical 
description of the slow amplitude dynamics, introduced 
an effective dynamical potential, incorporated the effects 
of thermal noise leading to tunneling between the po- 
tential wells, and extended our analysis to photothermal 
forces. Finally, we have pointed out that the intrinsic 
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hysteresis may be useful to implement a 'latching' scheme 
for detecting small displacements. The effects described 
here should be measurable in cavities that are within the 
reach of current technology. 
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